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Abstract 

We present a scheme in which an ion trap quantum computer can be used to 
make arbitrarily accurate measurements of the quadrature phase variables for 
the collective vibrational motion of the ion. The electronic states of the ion 
become the 'apparatus', and the method is based on regarding the 'apparatus' 
as a quantum computer register which can be prepared in appropriate states 
by running a Fourier transform algorithm on the data stored within it. The 
resolution of the measurement rises exponentially with the number of ions 
used. 
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Quantum computing offers promise for greatly enhanced efficiency of implementation for 
difficult computational problems The most important example to date is the quantum 

factoring algorithm of Shor |^ . Recent advances in quantum error correction and fault 
tolerant computation indicate that there is nothing in principle which will prevent the 
practical realisation of a quantum computer. The most promising suggestion currently is 
the ion trap realisation of Cirac and Zoller |T^. In this scheme, the internal electronic states 



of a string of ions becomes entangled with the lowest centre-of-mass vibrational mode of the 
ion trap. An experimental demonstration of a fundamental quantum logic gate in a trapped 



ion realisation was provided by Monroe et al [jTT 



In this paper we propose that an ion trap quantum computer be viewed as a means of 
making accurate measurements on the centre-of-mass mode of the ions. The electronic states 
of the ions forms an 'apparatus' coupled to the 'system', the vibrational mode. However us- 
ing the kind of discrete unitary transformations which realise a quantum computing circuit, 
we can prepare the apparatus in a variety of states to facilitate new quantum limited mea- 
surements on the system. This is in the spirit of the suggestion of Wineland and coworkers 
that highly entangled states may be used to improve frequency measurements. 

In recent experiments on trapped ions |T3] it was possible to subject the vibrational 



motion of the ion to a variety of unitary transformations by using sequences of Raman 
pulses. By selecting the Raman detuning carefully, a variety of oscillator states may be 
excited. In some case the unitary transformation of the vibrational motion is conditioned 
on the electronic state of the ion, and may be different for the ground and excited states. 
This conditional dependence was recently used to generate a Schrodinger cat state for the 
vibrational mode. In other cases the unitary transformation of the oscillator is independent 
of the electronic state. 

Consider a single ion laser cooled to the ground state of a ion trap in the Lamb-Dicke 
limit. Assume that the ion is illuminated by laser fields which are well detuned from reso- 
nance. Let the interaction Hamiltonian coupling the vibrational and electronic states of a 
single ion be given by 



Hi" ^ — h(aia^ + Q;*a)|e)i(e| (1) 

where o;, is proportional to a classical laser field acting on the ith-ion in a dispersive regime, 
a, are the creation and annihilation operators for the vibrational mode of the ion, and |e)j 
is the excited electronic state of the ith ion 

If we now generahse to the case of N ions in a hnear trap, the electronic state of the system 
may be described by a binary string S = {Si = l,0;i = 1, N) , assuming some ordering 
of the ions, and where a 1,0 represent the excited state and ground state respectively. 
Alternatively the state may be described by an integer k for which the binary string is the 
binary code k = J2iLi Si2^. There are 2^ possible states for the electronic system so that 
k — 0,1,...,K — 2^ — 1. Assume now that the coupling strengths ccj are adjusted from ion 
to ion so that ai — q;2\ The couphng Hamiltonian between the electronic states and the 
lowest collective vibrational mode of the ions may then be written 

Hi = h{aa} + a*a)T (2) 

where 

T^j:k\k){k\ (3) 

The phase of a can be adjusted, so for simplicity let us take it to be real. The unitary 
operator describing the coupling of the vibrational and electronic states at the end of a 
sequence of pulses is then given by 

where x — a + is the operator describing the in-phase component of the collective vibra- 
tional amplitude in an interaction picture rotating at the ion trap frequency, and r is a real 
parameter. 

It is clear that the interaction can realise a measurement of the vibrational quantity 
X. In this interpretation we regard the N ion register as an 'apparatus' to measure the 
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'system', the vibrational mode. However if all the ions are initially in the ground state, 
they will remain in the ground state under this unitary interaction as T commutes with 
the projection operators for each electronic energy eigenstate. In order for the interaction 
to change the state of the ions we must first pre-process the electronic register to a state 
which can be displaced by the unitary interaction. The required states \k) are eigenstates 
of an operator $ which is canonically conjugate to T. That is to say T must act as a pure 
differential operator in the {\k) k = 0, . . . , K} basis. This basis is simply a discrete Fourier 
transform of the original electronic basis {\k) k = 0, . . . , K}. Thus 

ir) 1 fexpfi?!^i)|.) (5) 



It is at this point that we remember the quantum computer interpretation of this system. 
The transformation of Eq (|^) may be realised by running a Fourier transform algorithm on 
the electronic register ions [p!^p!5 . 



Let the initial state of the system be 

l^) = l^).®|0) (6) 

where is an arbitrary vibrational state and |0) indicates all the ions are in the ground 
state. The first step is to apply a sequence of 7r/2— pulses to the ions to place them in a 
symmetric superposition of the ground and excited states. The resulting electronic state is 
precisely the state |0), the 0-state in the Fourier transform basis. This state is equivalently 
a uniform superposition over all possible electronic energy eigenstates. In the second step, 
the unitary interaction is then implemented to couple the vibrational and electronic states. 
In the third step an inverse Fourier transform is run on the electronic register. At the end 
of these three steps the state of the system is 

\ K rco 2-Kikl 

\m') = ——Y, dp4>{p)e- — \p + rk),^\l) (7) 

where is the momentum probability amplitude for the initial vibrational state. 

In the final step we readout the state of the electronic register. This is very much like 
the readout of the output register in the Shor algorithm. For trapped ions this can be done 



with very high quantum efficiency using quantum jump techniques [|T3|. The result of this 



readout is a binary string describing which ions are in the ground state and which are in 
the excited state. Equivalently the resuh is the integer I encoded by this binary string. The 
(unnormahsed) conditional vibrational state of the system after a readout result, I, is 

1^)'' = T^VT 2^ / dp(P{p)e ^^+1 \p + rk)^ (8) 
The probability for this result is found to be 

\P2/K , 2TTi(k-k')l 

Pil) = jf-Tyi E e-^^x{r{k - k')) (9) 



where 



X{k) = -= I dxe'^^P{x) (10) 



with P{x) being the initial distribution of the in-phase quadrature variable x, of the vibra- 
tional state. Clearly x(^) is just the characteristic function for the 'position' distribution 
for the initial vibrational state. In what follows we will assume that r = 1. 

The probability distribution for the output / is thus seen to a kind of discrete approxi- 
mation to the probability distribution P{x). In the case of a minimum uncertainty Gaussian 
state with position variance A, and zero mean, the readout distribution takes the form 

= TTTT + f ,<^' + ' - (ffr) ■ 

In the limit A ^ oo we have an effective uniform position distribution. In this limit it is 
easy to see that P{1) = 1 for all /. In the opposite limit A — we have a very well defined 
position at X = 0, and P{1) = 5ifi. In general the sum can be truncated at rritrunc = 
when subsequent terms contribute less than e~^, which we refer to as the tolerance e. For a 
large number of ions, the truncation becomes a good approximation. In figure |1| we illustrate 
the intermediate behaviour for = 9 ions and various values of A, (the tolerance used is 
e = 0.01. In these plots we have reflected the values at I > N/2 to negative values so as to 
ensure the plot is symmetric around 0. We see that as the number of ions increases we get 
a better and better approximation to the true 'position' distribution P{x). 



The dimensionless position x is related to the index Ihy x = 



I, thus the corresponding 



position distribution P'{x) obtained from P{1) is defined over the domain [0,27r). Thus for 
a fixed number of ions A^, the range of the uncertainty A in the initial position distribution 
P{x) is restricted, if the measured position distribution P'{x) is to be essentially the same 
as P{x). The lower limit imposed on A is determined by the number of ions in the register, 
whereas the upper limit on A is independent of the size of the register and is given by A^ax ~ 
10. As the uncertainty approaches this upper limit, the measured electronic distribution P{1) 
becomes fiat, so that it does not provide any more information about the actual position 
distribution P{x). However, we are assuming that the collective motion of the ions is in 
the Lamb-Dicke regime, and hence we require a priori that the uncertainty in the position 
is bounded by Ax <^ \fN /?], where 77 is the Lamb-Dicke parameter of the centre-of-mass 
vibrational mode. Alternatively, if we are given a lower limit for values of A, then we can 
determine the minimum number of ions Nmin required to be able to measure the initial 
position distribution P{x) accurately, and this is given by 



for a tolerance of e = 0.01 in the summation series i.e., this corresponds to Nmin = 25 for 
A = 10"^°, and N^^in = 10 ions for A = 0.1. 

The graphs plotted in Figures 0(a)- (b) demonstrate that the value of the uncertainty Ax 
in the measured position distribution P'{x) settles to the value of the initial uncertainty A 
as the number of ions increases past the appropriate minimum number Nmin, thus indicating 
that the measured position distribution P'(x) is approximately equal to P(x). 

We have illustrated a way in which an ion trap quantum computer can be used to mea- 
sure the distribution of the in-phase and out-of-phase quadrature variables for the collective 
vibrational motion of the ion. The method is based on regarding the ions as a quantum com- 
puter register which can be prepared in appropriate states by running a Fourier transform 
algorithm on the data stored in the ionic register. 

The first point to note is that we are not restricted to any particular quadrature phase 



A^, 



mm 




(12) 
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variable. As the phase of the laser pulses coupling the system and apparatus may be varied 
we can in principle measure any rotated quadrature phase. Given a sufficient set of such 
distributions a quantum tomographic data inversion could be done to reconstruct the initial 
vibrational state. 

The second point to note is that as the number of ions increases the states available 
to the computer rises exponentially, providing an increasingly accurate readout of a system 
variable with a continuous spectrum. This generalises previous schemes to readout vibra- 
tional states using only a single ion; a two state apparatus. In that case many repetitions of 
the measurement must be preformed to get an average over particular vibrational variables. 
In our scheme we can do better (albeit at the considerable expense of running a quantum 
algorithm), as we have an apparatus with many more states, and thus one which is better 
adapted to a vibrational degree of freedom. 

Finally we note that the scheme is capable of considerable generalisation. For example 
by choosing different couplings between the ion and vibrational modes we can get access 
to many other vibrational variables. For example if we choose a sequence of Raman pulses 
which provided a 'squeezing' interaction for the vibrational motion, we can measure the 
distribution of eigenstates of the squeeze operator W^. Furthermore, there may be other 



quantum algorithms that realise different measurement schemes. 
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FIGURES 

FIG. 1. A plot of the distribution for ion state readouts for an input Gaussian state in the 
centrc-of-mass coordinate, with two different variances A. The number of ions is 9. Dashed: 
A = 1.0. SoHd:A = 0.1 

FIG. 2. A plot of the estimated variance in the centre-of-mass coordinate versus ion number 
for two different values of the initial coordinate variance, (a) A = 0.1 (b)A = 1. 
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